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Abstract 

The post-post-Newtonian (2PN) accurate mass quadrupole moment, for com- 
pact binaries of arbitrary mass ratio, moving in general orbits is obtained 
by the multi-polar post Minkowskian approach of Blanchet, Damour, and 
Iyer (BDI). Using this, for binaries in general orbits, the 2PN contributions 
to the gravitational waveform, and the associated far-zone energy and angu- 
lar momentum fluxes are computed. For quasi-elliptic orbits, the energy and 
angular momentum fluxes are averaged over an orbital period, and employed 

to determine the 2PN corrections to the rate of decay of the orbital elements. 
PACS numbers: 04.25.Nx, 04.30.-w, 97.60. Jd, 97.60.Lf 
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I. INTRODUCTION 



Inspiralling compact binaries are one of the most promising sources of gravitational ra- 
diation for kilometer size laser interferometric gravitational wave detectors like LIGO |I[ 
and VIRGO [Q. The method of matched filtering will be employed to detect and extract 
information of binaries from the inspiral waveforms [0,0 • In this technique one cross cor- 
relates the noisy output of a detector with theoretical templates. For this technique to 
be successful, the templates must remain in phase with the exact - general relativistic - 
waveform as long as possible. If the signal and template lose phase with each other even 
by a cycle in the ten thousand as the waves sweep through the bandwidth of the detector 
their cross-correlation will be significantly reduced and one may lose the event altogether. 
Detailed works on data analysis aspects |5|-|8[ bears out this inference and one is forced to a 
description of the evolution of the binary system, using the best available theory of gravity 
to substantially higher accuracy than that provided by the lowest order Newtonian approxi- 
mation. The construction of accurate theoretical templates for inspiralling compact binaries 
involves the solution to two different but related aspects referred to respectively as the "wave 
generation problem" and the "radiation reaction problem". In the generation problem one 
computes the gravitational waveforms and the associated energy and angular momentum 
fluxes emitted by the binary for a fixed, specified orbital motion ignoring the back reaction 
of the radiation emission on the orbit. In the radiation reaction problem on the other hand, 
one computes the effect of the emitted radiation on the orbital phase evolution and this is 
of crucial importance as explained earlier. 

Einstein's |§ far field quadrupole equation is the solution to the generation problem to 
the lowest order but applies only to objects held together by non-gravitational forces. Fock 



PCI and Landau-Lifshiftz [l lj provided two very different methods to generalize the above 
results to weakly self-gravitating systems and the two approaches are the starting points 
for the two methods available today to calculate gravitational wave generation to higher 
orders: the Blanchet Damour Iyer(BDI) ]T2| approach and the Epstein Wagoner Thorne 
Will Wiseman(EWTWW) [0 approach. 

The BDI approach builds on a Fock type derivation using the double-expansion method 
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of Bonnor fl4|| . This approach makes a clean separation of the near- zone and the wave- 
zone effects. It is mathematically well-defined, algorithmic and provides corrections to the 
quadrupolar formalism in the form of compact support integrals or more generally well- 
defined analytically continued integrals. The scheme has a modular structure : the final 
results are obtained by combining an 'external zone module' with a 'near zone module' 
and a 'radiative zone module'. For dealing with strongly self-gravitating material sources 
like neutron stars or black holes one needs to use a 'compact body module' supplemented 
by an 'equation of motion module' to describe their 'conservative' orbital motion. Using 
this approach the generation of gravitational waves from inspiralling compact binaries of 
arbitrary mass ratio moving in a quasi circular orbit has been computed to 2PN accuracy 



p5| . [T6|| and more recently to 2.5PN accuracy ||17|| . In this paper, in the first instance, using 
the BDI approach, we extend the above 2PN treatment to inspiralling compact binaries 
moving in a general orbit and compute the 2PN contributions to the waveform and the 
energy flux. Unlike for circular orbits, the angular momentum flux from general orbits 
provides additional, independent information and we also compute the same. 

The Epstein- Wagoner-Thorne- Will- Wiseman approach on the other hand, builds on a 
Landau-Lifshiftz type treatment to derive post- Newtonian corrections to the lowest order 
quadrupole formula. The combined use of an effective stress energy tensor for the gravita- 
tional field ( with non-compact support ) and of formal post-Newtonian expansions led to 
the appearance of divergent integrals. The presense of the divergent integrals and the lack 
of a clear separation between the near zone and the wave zone were unsatisfactory features 
of this scheme until recently. However, last year, Will and Wiseman |18| have provided a 
resolution to this problem by taking literally the statement that the solution is a retarded 
integral i.e., an integral over the entire past null cone of the field point. A careful evaluation 
of the far-zone contributions, then shows that all integrations are indeed convergent and 
finite and moreover the tail terms are also correctly recovered. Using this treatment, Will 
and Wiseman have computed the 2PN accurate waveform and energy flux for general orbits. 
We thus have two approaches to the 2PN generation, which can provide a useful check on 
the long and tedious algebra. 

The most accurate results to date for the generation and the radiation reaction have been 
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obtained in the limit where a test body orbits a very massive central body. In this comple- 
mentary approach, based on black hole perturbation techniques there exist numerical results 
that are exact in (v/c) and analytical results accurate to the 5.5PN order, i.e. corrections 
of 0[(v/c) u ] for a test particle in a circular orbit around a Schwarzschild black hole 



where v is the orbital velocity of the test particle. For a test particle in a slightly eccentric 
orbit, around a spinning black hole expressions for the energy and the angular momentum 
fluxes have also been computed to the 2.5PN order [24|. 

It is well-known that, when gravitational waves from prototype systems like the binary 
pulsar 1913+16 enter the bandwidth of the terrestial interferometric detectors, the eccentric- 
ity of these binary systems would have been drastically reduced and have become negligible 
due to the gravitational radiation reaction. A treatment of such systems is simpler since 
the quasi-circular approximation for their orbits is amply adequate and the corresponding 
waveforms do not depend on the eccentricity parameter of the orbit. However there exist 
scenarios in which the eccentricity is no longer negligible and this would require the more 



general treatment provided in this paper and independently by Will and Wiseman |Tj| . One 
such possibility has been discussed by Shapiro and Teukolsky in the context of the for- 
mation of supermassive black holes. They consider a cluster of compact objects -neutron 
stars and black holes- residing at the center of a galactic nucleus. Coulomb scattering and 
dissipative processes will drive such a cluster to a high density, high redshift state. Once 
the central redshift becomes sufficiently large, relativistic instability sets in, and the core 
undergoes catastrophic collapse to form a supermassive black hole. Quinlan and Shapiro 
p8| have shown that during the final year of the evolution of such a cluster, just prior to the 
catastrophic collapse, there can be 100 — 10 4 evolving black hole binaries in eccentric orbits 
driven by gravitational radiation reaction, with masses in the range 10 — 100 M Q . Ground 
based interferometric detectors will be sensitive to the gravitational radiation from these 
binaries in eccentric orbits and such eccentric binaries may prove to be another possible 



class of gravitational wave sources. More recently, Flanagan and Hughes |29| suggest that 
intermediate mass black hole binaries of the kind considered above -with total masses in 
the range 50M Q < M < (a few) x 1O 3 M - may well be the first sources to be detected by 
LIGO and VIRGO. The other possibility involves compact objects orbiting 10 6 to 1O 7 M 
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black holes, that seem fairly common in galactic nuclei. In this case the compact objects 
could be scattered into very eccentric orbits orbits via gravitational deflections by other 
stars. However, by the time gravitational radiation reaction becomes the dominant orbital 
driving force, there is not enough inspiral remaining to fully circularise these orbits. Hils and 
Bender [^UJ have argued that the event rates for the above process are very encouraging and 
the chances of such signals being observed by Laser Interferometric Space Antenna, LISA 



3jJ appear very good. 

The expressions for the far-zone energy and angular momentum fluxes find application 
in another related but distinct problem; the evolution of the orbital elements of systems 
like the binary pulsar 1913+16; most importantly the orbital period, and to a lesser extent 
the eccentricity and semi-major axis. This application, in addition to the generation results 
discussed earlier, requires a convenient representation of the post-Newtonian motion of two 
point masses in elliptical orbits. To 1PN accuracy, such a quasi-Keplerian representation 
has been provided by Damour and Deruelle |32j , while to the 2PN order a generalised quasi- 



Keplerian representation has been implemented by Damour, Schafer, and Wex [|3~3| -|3~6"|j . This 
representation differs from the Keplerian representation of the Newtonian motion through 
the appearance of three eccentricities instead of one, and a constant measuring the secular 
advance of the periastron. Starting from the above representation of the orbital elements 
in terms of the conserved energy and angular momentum, one computes the time variation 
of the orbital elements. One ends up with a result, in terms of the time variation of the 
'conserved' energy and angular momentum. By a heuristic argument, one replaces these by 
the corresponding average far- zone fluxes which may be computed by averaging the far-zone 
fluxes over an orbital period, using the quasi-Keplerian orbital representation. The reduction 
in the orbital period, accurately inferred from the timing data of the binary pulsars is in 
excellent agreement with the rigorous predictions of general relativity [p7]-|4l|, which in turn 
are consistent with the results of the above heuristic approach |42j44|,^5|. Extending the 
above approach, Blanchet and Schafer have obtained the 1PN and the 1.5PN corrections 
to P, the rate of decay of the orbital period P [f5],[17]]. They have shown that for PSR 
1913+16, the relative 1PN and 1.5PN corrections are numerically equal to +2.15 x 10 -5 
and +1.65 x 10 -7 respectively. These are unfortunately far below the present accuracy in the 
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measurements of P for 1913+16. Junker and Schafer [|8| computed the 1PN contributions 



to the gravitational waveforms, the associated angular and linear momentum fluxes and used 
it to compute the evolution of the orbital elements in the quasi-Keplerian representation. In 
the other part of the paper, we extend the above computations to obtain the 2PN corrections 
to the evolution of orbital elements, taking due care of a new complication at this order that 
the far- zone fluxes are computed in the harmonic or De Donder coordinates, while the orbital 
representation is available in the Arnowit, Deser and Misner (ADM) coordinates. 

Briefly, in this paper we obtain the terms 0(e 2 ) in the expressions below, where e ~ 
v 2 /c 2 ~ Gm/c 2 r; m,r,v, being the total mass, the distance between the bodies and the 
relative velocity of the two bodies, respectively, 



and where iy is the mass quadrupole moment for a system of two compact objects moving 
in general orbits while is the transverse-traceless(TT) part of the radiation-field, repre- 
senting the deviation of the metric from the flat spacetime. In the above, d£/dt , dj7"/dt are 
the far- zone energy and angular momentum fluxes, < dS/dt > and < dj jdt > represent the 
averages of the far-zone fluxes over an orbital period, while < da r /dt >, < de r /dt > along 
with dP/dt give the gravitational radiation driven rate of decay of the orbital elements of 
the binary in the generalized quasi- Keplerian parameterization. Note that the suffix 'N' 
denotes Newtonian contribution in all the above equations. For example (/^tDn denotes the 




(1.1a) 
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Newtonian contribution to the waveform given by {2 G / (c 4 R)} Vijkm ifj ■ See Eq. ( |5.1|) for 
our notation. 

The plan of the paper is as follows: In section [TT], using the BDI approach, we compute 
the 2PN accurate mass quadrupole moment for two masses moving on general orbits. We 
also obtain and list the mass and the current moments to the required accuracy, needed to 
compute the 2PN accurate waveform. In section III we calculate the 2PN contributions to 
the far-zone energy and angular momentum fluxes and discuss the limiting forms of these 



expressions. In section [TV], for the quasi-elliptic case, we average the above expressions over 
an orbital period to obtain the 2PN corrections to < dS/dt >, < dj jdt > and the rate of 
decay of the orbital elements. Section [V| computes the 2PN contribution to gravitational 
waveform for general orbits. Section |VI] contains the summary and a few concluding remarks. 
In Appendix A we list identities, that are used in the computations, especially, of the 
waveform. Finally, in appendix B, we sketch the steps involved in verifying the equivalence 
of our waveform obtained using the STF multipoles of the radiative field and the Will- 
Wiseman one obtained using the Epstein Wagoner multipoles. 



II. MASS AND CURRENT MOMENTS OF COMPACT BINARIES ON 
GENERAL ORBITS FOR 2PN GENERATION 

A. 2PN mass quadrupole moment 



The starting point for the computation of the 2PN accurate mass moment is the form of 
the moment quoted in []I6| i.e., Eq.(2.17). 
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The symbol FPb=o m the above stands for "Finite Part at B = 0" and denotes a mathe- 



matically well-defined operation of analytic continuation. For more details see [|T6 [. 

As emphasized in though the above expression is mathematically well defined, it is a 
non-trivial and long calculation to rewrite it explicitly in terms of the source variables only. 
This is achieved by representing the stress energy tensor of the source as a sum of Dirac 
S-f unctions. 

T ^ t) = £ m /ff^p {x _ yA(t)h (2 . 2) 

where m A denotes the (constant) Schwarzschild mass of the Ath compact body. Evaluating 
this to 2PN accuracy we obtain for the source variables 

/ v 2 



a(x, t) = £ fi A (t) 1 + -f U(x - y A (t)) , (2.3a) 

A=l \ C / 

N 

vifct) = ^MW-yAW) , (2.3b) 

A=l 

N 

cr y (x,t) = J2 ^Ait)v A v A 5(-K-y A (t)) , (2.3c) 

A=l 

where v\ = dy A /dt and 

fi A {t) = m A {l + (d 2 ) A + (d^A , (2.4a) 



1 n , 



rf 4 = I |^v 4 + ^f/v 2 - 4Cfo - 2$ + ^f/ 2 + 4t/ ss } , (2.4c) 

In the above V denotes the combination 

V = U+±- 2 d 2 t X, (2.5) 

the potential appearing naturally in the 1PN near- zone metric in harmonic coordinates. The 
subscript A appearing in Eq. ( |2.4a| ) indicates that one must replace the field point x by the 
position y A of the Ath mass point, while discarding all the ill-defined (formally infinite) 
terms arising in the limit x — > y A . For instance 

(t/^GE ^f'y . (2.6a) 

B+A IYA-YbI 
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[Note that the second time derivative appearing in V, Eq. (2J3), must be explicated before 
making the replacement x — > yA{t).} 

The terms in Eq.( [2.1|) fall into 3 types: compact terms, Y terms and W terms. The 
compact terms, where the 3-dimensional integral extends only over the compact support of 
the material sources; the Y terms involving three dimensional integral of the product of two 
Newtonian like potentials; and the W term involving three dimensional integrals of terms 
trilinear in source variables. The evaluation of these different terms proceeds exactly as in 
the circular case. In fact if the time derivatives are not explicitly implemented the expression 
in the general case and the circular case would be identical. The difference obtains when the 
time derivatives are implemented using the equation of motion. In this section we need to 
use the general form of the Damour-Deruelle equations of motion rather than the restricted 
form of the circular orbit equations of motion relevant in [|16[ . 



We take up the compact terms first. They are given by 
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in which we have introduced for convenience jl A = //a(1 + v a/ c2 )- m ^ e above form the 
moment depends not only on the position and velocity of the bodies but on higher time 
derivatives. It is in the reduction of these derivatives that we need the 2PN accurate equation 
of motion for general orbits. We use a harmonic coordinate system in which the 2PN center 
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of mass is at rest at the origin. Using the 2PN accurate center of mass theorem, in the 
center of mass frame, we can express the individual positions of the two bodies moving in 
general orbits in terms of their relative position x = y x — y 2 and velocity v = v x — v 2 
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where r = |yi — y 2 | is the harmonic separation between the 2 bodies. The explicit values 
of Xi an d X2 are not needed in our calculations and hence not given above. The above 
equations are obtained by setting equal to zero the conserved mass dipole G for general 
orbits. Here we denote 



m = mi + m2 , 5m = m\ — m 2 , 

X lS ?L, X 2 = ™ =1 _ Jfl , 

m m 
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m^ m 

The 2PN accurate equations of motion is written down next for completeness, where 
finite-size effects, such as spin-orbit, spin-spin, or tidal interactions are ignored. |49|-|5T|. For 
the relative motion we have 



a = a N + a£L + al^L + 0(e 5 ) 



(2.10) 



where the subscripts denote the nature of the term, Newtonian (N), post-Newtonian (PN), 
post-post-Newtonian (2PN), and the superscripts denote the order in e. The explicit ex- 
pressions for various terms mentioned above are given by 
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where n = x/r and f = dr/dt. 

We have on hand all the ingredients to compute I c . Though long and tedious the com- 
putation is straightforward and yields for the 2PN mass quadrupole: 
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The Y terms on the other hand are given by 
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The explication of all the above terms finally leads us to 
ry] 2mr] Gm 
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The evaluation of the 1^ term, the new feature at 2PN level, was discussed in detail in 
||16|| . The W term has been evaluated there for general orbits and we need to use the same 
result here. We have 
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Adding up the compact i.e., C, Y and W contributions given by Eqs. (|2.12|) , (|2.16|) and 



( 2.17 ), we finally obtain the expression for the 2PN accurate mass quadrupole for a system 
of two bodies moving in general orbits. The final result is written below as a combination 
of the three possible combinations x lJ ,x l v J ,v lJ with coefficients which include corrections 
beyond the Newtonian order at 1PN and 2PN orders: 
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(2.18) 



The above expression is identical to the one obtained by Will and Wiseman in the appendix E 
of using the new improved version of the Epstein- Wagoner formalism. In their treatment 
the Epstein- Wagoner multipoles appear more naturally, using which they compute the STF 
mass quadrupole moment. Since the approach employed here and in jTj| follow algebraically 
different routes, the above match provides a valuable check on the long and complicated 
algebra involved in the determination of the crucial mass quadrupole moment for 2PN 
generation. 



B. The other relevant mass and current moments 



In this section we list the higher order mass and current multipole moments, required to 
compute the 2PN contributions to the gravitational waveform and the associated far-zone 
energy and angular momentum fluxes. They are straightforwardly obtained by explicating 
the point particle limits of the more general expressions in the earlier BDI papers JT2| 
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The mass and the current moments listed above, agree with Eqs. (E3) of For the case 



of circular orbits, the above mass and current moments reduce to Eqs. (4.4) of [jnj 
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III. THE FAR-ZONE FLUXES 



A. The Energy flux 



As discussed in the end result of the 2PN accurate generation formalism is an 
expression relating the radiative mass and current multipole moments U l and Vl respectively 
to the source mass and current multipole moments II and Jl respectively, obtained in the 
previous section. In particular, 



U l3 {T R ) 
Uijk(T R ) 
Vij(T R ) 



2Gm 



dr 



C JO 

2Gm r+o. 



dr 



C° JO 

2Gm r+o. 



dr 



ln S) + i 

\2b) 6 



lU(T R ~ r) + 0(e 5 



4?(T R -T) + 0(e*) 



C° JO 

for the moments that need to be known beyond the 1PN accuracy, and 



I®(T R )+0(e*) , 



Ul(T r ) - * L 
V L (T R ) = jP(T R )+0(e 3 



(3.1a) 
(3.1b) 
(3.1c) 

(3.2a) 
(3.2b) 



for the other ones. The integrals in the above expressions, associated with the gravitational 
wave tails, contain in addition to the total mass-energy of the source m, a quantity b which 
is an arbitrary constant with dimensions of time parametrizing a freedom associated in the 
construction of the far-zone radiative coordinate system. 

In terms of the STF radiative moments of the gravitational field the far- zone energy flux 
to 2PN accuracy is given by |7| (with = d n U/dT%): 



'd£T 
dt 



^/V^ (1) + - 



far— zone 
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l 



— ^— /7 (1) /7 (1) 4- — V-^'V- 

1 on ^ijk^ijk i AK V l 3 



1 



189 
1 



45 



9072 % i km ^ km 84 l i k ^ k 



r(l) T /(l) 



Oie" 



(3.3) 



The 2PN-accurate energy loss given by Eq. (|3.3|) splits naturally into an "instantaneous" 
contribution and a "tail" one. In this paper, we deal only with the instantaneous contribu- 
tion, which is given by [0,0 

■ 1 



dt ) f 

\ / tar— zone 



^/i/( 3 )/.( 3 ) 



: 5 1 5 



16 



(3) t(3) 

45 _ y v 



+ 



1 r(5) r(5) 

9072 j km ijk m 



1 

84 



t(4) t(4) 
" ijk' J ijk 



(3.4) 
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Here denotes the n th time derivative of STF multipole moment of rank L. All the 



computations from here onwards are performed, using MAPLE |53| . Evaluating the relevant 
time derivatives of the multipole moments in Eq. fl3.4| ), using the post-Newtonian equations 
of motion to the appropriate order we obtain 



-jT) -t N + t 

III / far— zone 



1PN + &2PN 

8 G :i m 2 fi 



£ 



1PN 



{l2v 2 - llr 2 } , 

(785 - 852t/)t; 4 



15 c 5 r 4 
8 G 3 m V f 1 
15 c 7 r 4 1 28 

- 2(1487- 13927/)t> 2 r 2 

- 160(17 — 77) v 2 

r 

+ 3(687 - 620r/)r 4 + 8(367 - Ih-q) 



Gm 



+ 16(1 -4?7) 



G 2 m 2 



2PN 



8 G 3 m 2 /i 2 

15 c 9 r 4 
1 



42 



;i692 - 5497r/ + 4430^)u 



- —(1719 - 10278r/ + 6292t/ 2 )t; 4 7 ; - 2 

1 Gm 

(4446 - 52377/ + 1393t/ 2 ) tj 4 

21 ' r 

+ ^(20 1 8 - 152 77/ + 75 727/ 2 )7; 2 f 4 

1 Gm 

+ -(4987 - 85137/ + 2165t/ 2 ) i; 2 f 2 

7 r 



+ —(281473 + 81828// + 4368r/ 2 
756 



G 2 m 2 



(2501 - 202347/ + 8A0Ari 2 )r 6 



1 

42 

— (33510 - 609717/ + 14290t/ 2 ' 
63 

— (106319 + 97987/ + 5376t/ 2 ; 
232 



r 

G 2 m 2 



+ —(-253 + 10267/ - 56r/ 2 
63 



G 3 



m 



(3.5a) 
(3.5b) 



(3.5c) 



(3.5d) 



Eqs. Q3.5| ) are in exact agreement with the results of Will and Wiseman using the new 
improved Epstein- Wagoner approach |TH] . Circular and radial infall limits of Eqs. (|3.5|) are 



in agreement with earlier results 



ctO| , |T8"[ and discussed further in section |111 C| . 
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The tail contribution, on the other hand is given by 

tail 



'd£_\ 
dt) 



2G 2Gm 



5c 5 



R-T) 



(3.6) 



far— zone 

where b\ = 6e~ n / 12 . A detailed discussion of the tail terms and its implications has been 



given by Blanchet and Schafer [47], and we do not discuss it any further in this paper. 



B. The angular momentum flux 



In terms of the STF radiative multipole moments the far-zone angular momentum flux 
to 2PN accuracy is given by: 

G 



dj 
dt 



far— zone 



2 U .[/W + I 



1 TT ttO) i 32 V T/« 



1 



2268 



1 

28 



(i) 



(3.7) 



As before, rewriting the radiative moments in terms of the source moments, allows us to 
separate the instantaneous and tail contributions and we discuss them independently. We 
have IE21 



dj 
dt 



inst 



far— zone 



/^jO) + I 
1 



J_ r(3) r(4) 32 (2) (3) 
gg Pjk A qjk 45 PJ 



' r(4) r(5) no; 
2268 3 kl-Lqjkl 2g Pjk J qjk 



r(3) t(4) 



(3.8) 



Computing the required time derivatives of the STF moments, using the post-Newtonian 
equations of motion to the appropriate order, we obtain 



V dt / far- 



Jn 

JlPN 



Jn + JlPN + ^2PJV 

8G 2 m/i 2 

5 c 5 r 3 
8G 2 m/i 2 



L N { 2v 2 - 3r 2 + 2 



Gm 



(3.9a) 
(3.9b) 



5 c 7 r 3 " Ln 184 



(307 - 548r))v 4 
Gm 



- 6(74 - 277r7)t> 2 r 2 - 4(58 + 95rj) 
+ 3(95 - 360r/)r 4 + 2(372 + 197/7) 



r 
Gm 

r 



2(745 - 2-q) 



G 2 m 2 ' 



(3.9c) 
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• _8G 2 m// 2 ~ 

U2PN — n — 



1 



(2665 - 123557/ + 1289 4t/ 2 )tj 6 



5 c 9 r 3 " 1504 
3(2246 - 1265377 + 156377/ 2 )?/f 2 



+ (165 - 491/7 + 4022t/ 2 



Gm 



+ 3(3575 - 168057/ + 15680r] 2 )v 2 f 4 



+ (21853 - 216037/ + 255177 2 



Gin 



, v r 



G 2 m 2 

- 2(10651 - 101797/ + 3428T7 2 ) — — v 2 



28(195 - 815t/ + 4857/ 2 )r 6 



(22312 - 41 398 7/ + 96 95 t/ 2 )?^ 4 



Gm 



G 2 m 2 

+ 2(8436 - 251027/ + 4587t/ 2 ) — — r 2 



+ 



2268 



(170362 + 704617/ + 1386t/ 5 



G 3 m 3 



(3.9d) 



where Ln = r x v. The 1PN contribution is in agreement with the earlier results of Junker 
and Schafer [IS ] . The 2PN contribution is new and together with the energy flux obtained in 
the earlier section forms the starting point for the computation the 2PN radiation reaction 
for compact binary systems -4.5PN terms in the equation of motion - 



using the refined 



balance method proposed by Iyer and Will 5^.57]. The tail terms, in the angular momentum 
flux are given by 

tail 

p R -r). (3.10) 



'dJi 



dt 



2G 2Gm (3) 

€ ijk hi 



5c 5 c 3 



far— zone 

We will not be discussing the tails terms here, as they are extensively studied by Rieth and 
Schafer 



C. Limits 

All the complicated formulae, discussed in the earlier sections take more simpler forms for 
quasi-circular orbits. For compact binaries like PSR 1913+16, quasi-circular orbits should 
provide a good description close to the inspiral phase, since gravitational radiation reaction 
would have reduced the present eccentricity, to vanishingly small values. In this context 
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'quasi' implies the slow inspiral caused by the radiation reaction. The quasicircular orbit is 
characterised by f = f = 0(e 2 ' 5 ). The 2PN equations of motion, become 

dv g? 2 x 



dt dt 2 

with co>2PN; the 2PN accurate orbital frequency, is given by 
o Gm ( ( 41 



u4 N x + 0(e 2 - 5 ) . (3.11a) 



10 



2PN - ~ 



{l - (3 - 77)7 + (6 + + r? 2 ) 7 2 | , (3.12) 



where 7 = Gm/c 2 r. Note that Eqs. (|3.11|) imply as usual, that v = |v| = W2PN r + 0(e 2 - 5 ) 
so that from Eq. ( |3.12j ) we get 



v 2 = ^ {l - (3 - 77)7 + (6 + ^ + r/ 2 ) 7 2 } • (3.13) 

Substituting r = in Eqs. (|3.5|), and using Eq. (|3.13|) we obtain the 2PN corrections to the 
far- zone energy flux for compact binaries of arbitrary mass ratio, moving in a quasicircular 
orbit 

fd£\ inst 32c 5 , . r /2927 5 \ , / 293383 380 \l /0 _ x 

U , = ^g" 7 T-nw+f +7 boTT+if ■ (314) 

\ / tar— zone 

Eq. ( |3.14| ) is consistent with results of 



The energy and angular momentum fluxes are not independent but related in the case 
of circular orbits. The precise relation may be written following as: 

{'dt) 

where (^f-) = L Nl j , (3.15b) 

V dt J far— zone 



= v 2 J, (3.15a) 

far— zone 



where v 2 defined in terms of Gm/r, is given by Eq. (|3.13|) [16 



The other limiting case we compare to corresponds to the case of radial infall of two 
compact objects of comparable masses. Equations representing the head-on infall can be 
obtained from the expressions for the general orbit by imposing the restrictions, x = zh, 
v = zh, r = z and v = r = z. We consider two different cases, following Simone, Poisson, 
and Will [pDJ. In case (A), the radial infall proceeds from rest at infinite initial separation, 
which implies that the conserved energy E(z) = E(oo) = 0. In case (B), the radial infall 
proceeds from rest at finite initial separation zq, which implies 
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E(z = zq) = -fic 7o i 1 - o 7o + o 



1 



1 



i 15 

H T) 

2 ' 



7o 2 



2 ' 2 

where 70 = Gm/ ZqC 2 . Inverting -E(z) for i: 2 and using Eq. ( |3.1£ ) we obtain 



(3.16) 



z = -c \ 2 (7 - 70 ) 



l--(l-|)+7o(l-^ 



+7 2 (l3-^ + 5r ? 2 



4 



"„(r,-^ + lV 



57/ 



(3.17) 



where 7 = Gm/zc 2 . Using the radial infall restrictions and Eq. ([3 -IT ) in Eqs. ( |3.5|) we 
obtain for the case B),the far-zone radiative energy flux 

inst 



~dt 



far— zone 



16 
15 



c 5 rf 7 5 U 



,r 



111 / \ 9/ 135 7? , 
43 - —77 - x(116 - 13177) + x 2 (71 — 



7 



1127 803 2 x /4471 15481r? A 2 

r/- 112r/ 2 + - + 2864?7 2 



9 



12 



7 V 9 



x" I 38521?7 t 8800/7^ \ 3 1183/7 , 872^ 



— I 1870 



+ 



+ x 6 83 



+ 



7 



6 3 J \ 4 7 

where x = 70/7. For the case (A), the expressions for z and dS/dt are obtained by setting 
7o = in Eqs. ( p.!7|) and ( |3.18| ). Eq. ( |3.18j ), along with corresponding one for case (A) 
are in agreement with |60|j. 



(3.18) 



IV. THE EVOLUTION OF THE ORBITAL ELEMENTS IN THE GENERALIZED 
QUASI-KEPLERIAN PARAMETERIZATION OF THE BINARY 



In this section, we compute the 2PN corrections to the rate of decay of the orbital 
elements of a compact binary, in quasi-elliptical orbit, i.e the effect of the 4.5PN radiation 
reaction on a 2PN accurate conservative elliptical motion, extending the earlier computations 
4"6|-4"8" . The basic ingredients we employ for the calculations are the far-zone energy and 



angular momentum fluxes in the harmonic coordinates computed in previous sections and 
a 2PN accurate description of the relative motion of the compact binaries available in a 
generalized quasi-Keplerian parameterization given in the ADM coordinates P^-P^j. Since 
the De Donder(harmonic) and the ADM coordinates are different at the 2PN order, we 



use the coordinate transformations connecting the harmonic and the ADM coordinates [61 
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to rewrite the far-zone fluxes in the ADM coordinates. The far zone fluxes, in the ADM 
coordinates are averaged over an orbital period, extending the earlier computations at the 
1PN and the 1.5PN order [p]|- f4l||5B|j . The 2PN corrections to the rate of decay of the orbital 
elements are computed using heuristic arguments based on the conservation of energy and 
angular momentum to the 2PN order. Before proceeding to the actual computations, in the 
next two sections, we summarize the generalized quasi-Keplerian description of the binary 
orbits in the ADM coordinates and the transformations needed, to relate the kinematical 
variables in the harmonic and the ADM coordinates. 



A. The second post- Newtonian motion of compact binaries 

The generalized quasi-Keplerian description for the general binary orbits to the 2PN or- 
der, developed by Damour, Schafer, and Wex p4"H36[1 is best suited for the calculation we pro- 



pose to do in the following sections and we summarize it in what follows. Let r^(t^), c^a^a) 
be the planar relative motion of the two compact objects in usual polar coordinates asso- 
ciated with the ADM coordinates. The radial motion ta^a) is conveniently parameterized 
by 

r A — a r (l — e r cosw) , (4.1a) 

n(tA — to) = u — tt sin u -\- —- sin v -\- —j (v — u) , (4.1b) 

c 4 c 4 

where is the 'eccentric anomaly' parameterizing the motion and the constants 

and to are some 2PN semi-major axis, radial eccentricity, time eccentricity, mean motion, 

and initial instant respectively. The angular motion 4>A(tA) is given by 

0a-0o= U + ^U + ^sin2u + ^sin3u, (4.2a) 



where v = 2 tan" 1 { — tanQ } . (4.2b) 




In the above <fto, k, e$ are some constant, periastron precession constant, and angular ec- 
centricity respectively. All the parameters n, k, a r , e t , e r , e^, f t , gt, and are functions 
of the 2PN conserved energy and angular momentum per unit reduced mass/z: To avoid 
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additional notations following fl34|-|36]j , these are also denoted as E and |J|. Their explicit 



functional forms, given in [55 are displayed below 



Gm 

Ye 



i(l + 10^ + ^)£ 2 + i(17-ll,,)^ 



n 



1 + 2Eh 2 - ^{2(6 - rj)E + 5(3 - 77) E 2 h 2 ^ + ^{(26 + 7/ + rf)E 

E ~\ 
-2(17 - II77) — + (80 - 55t/ + 4r] 2 )E 3 h 2 \ , 

{-2E) 



Gm 



4 c 



' (15-^+1 



1(555 + 307/ + 11 V )E 2 - ^(5 - 277) { -^- 



1 + 2Eh 2 + 1|4(1 - rj)E + (17 - 7 V )E 2 /i 2 } + ^2(2 + 77 + 5 V 2 )E 2 



-(17 — II17) — + (112 - 477/ + 16r] 2 )E 3 h 2 - 3(5 - 2t/)(1 + 2£/Y 



-2£) 



h 2 



ft = -^j i ^ + r l ){l + 2Eh 2 Y( 
r/ t =-(5-27/)- 1 



-2E) 



h 



h 2 
1 7/ 

8 



2 c 2 



(5-2t/)£ + — (7-2r/) 



(1 -3t/)(1 + 2E/i 2 



32 /i 4 v ; ' 

1 + 2Eh 2 - ^ll2E + (15 - T])E 2 h 2 \ - — r (4(16 - 8877 - 9t/ 2 )£ 2 
<r I J 8 c 4 I 

-4(160 - 30t/ + 3r] 2 )E 3 h 2 + (408 - 232// - 15r/ 2 ) — ^ , 

/V J 



(4.3a) 

(4.3b) 
(4.3c) 

(4.3d) 
(4.3e) 

(4.3f) 
(4.3g) 
(4.3h) 
(4.3i) 



(4-3j) 



where h = |J|/(Gm). Note that n = 2tt/P, where P is the period of the binary. Using 
these parametric equations of the motion, we compute r\, v\ to the 2PN order in terms of 
E, h 2 , (1 — e r cos u) using, 



du 


dt^ dtA dv 
du dv du 


(4.4a) 


r\ = 


/ dr A ^dt A \ 2 
\ du du J 


(4.4b) 


& = 


f d(j)A dv ! dtA \ 2 
y dv du du J 


(4.4c) 




r A + r 2 A <p 2 A- 


(4.4d) 
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The subscript 'A ' present in Eqs. ( |4.4j ) is a reminder that the expressions refer to the ADM 
gauge. We have 

2 2 



-1 + 



+ 



(1 — e r cosw) (1 — e r cosii) 2 



Eh 2 (-2E) 



+ ^(-3 + 9^ + 
1 

(1 — e r cosu) 



1 — e r cosu) 



38 - 30r/ 



— - 4 - 1977 + 16rf 



40 - 20t/ - (36 - 28t/) Eh 2 
1 



1 



1 — e r cosu) 



(64 - 24t/) Eh/ 



E 2 



(1 — e r cosu) 



168 - 326r/ + 98t/ 2 - -— (34 - 22t/) 



'1 — e r cosu) 2 
1 



496 - 712t/ + 164t/ 2 - (213 - 298t/ + 85t/ 2 ) £/i 2 
212 - 332t/ + 80t/ 2 - (800 - 932t/ + 188r/ 2 ) £/i 2 



1 — e r cost/) 3 

' 528 - 528t/ + 96t/ 2 



v l — e r costj 



£/i 2 + 



1 



1 — e r cosu) 



32 + 8t/ 2 t/£ 2 /i 4 }>(-£) 3 , 
(4.5a) 



1 — e r cosu) 



}(-2E) 

- 97/ - [38 - 30t/] + 

c z I (1 — e r c~ ■ ■ 



r COS U ) 
2 I v2 



1 — e r cosu) 2 



[40 - 20t/] 



+8 71 " V^ Eh2 ) E 

(1 — e r cosu) 6 J 



— 4- 19t/ + 16t/ 2 



(1 — e r cost/) 



168 - 326t/ + 98t/ 2 



Eh 2 



(34 - 22t/) 



'(1 


— e r cosu) 2 
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(1 


— e r cosn) 3 




1 



428 - 6687/ + 164t/ 2 

212 - 332t/ + 80t/ 2 - (76 - 84t/) t/ E/i 2 

80 - 128t/] 7/ £/i 2 + 72 rf E 2 h 4 \(-Ef 

(1 — e r cosur ) 



(1 — e r cosm) 4 

These expressions for r\ and v\ are consistent with Eqs. (6) and (7) of 



(4.5b) 



B. The transformation between De-Donder ( harmonic ) and ADM gauges 

As pointed out earlier, the far-zone fluxes obtained in previous sections are in the har- 
monic coordinates, whereas, the 2PN accurate orbital description given by Eqs. ( |4.1|) , ( (4.2|) , 
and ( |1.3|) are in the ADM coordinates. For the purpose of averaging the far zone fluxes using 
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the the 2PN accurate orbital representation, we need to go from the De-Donder (harmonic) 
to the ADM gauge, and rewrite the expressions for the far-zone fluxes in the ADM coor- 
dinates. These follow straightforwardly from the transformation equations in [EJI| and we 
list below the transformation equations, relating the harmonic(De-Donder) variables to the 
corresponding ADM variables: 



r D = r A + 



Gm 
8c 4 r 



(5v 2 - r 2 ) 77 + 2 (1 + 127?) 



Gm 



tv — t\ 



18 77 rrv} , 

Gm 
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8c 4 r 2 
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c 4 r 



(Lat)d 
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Gm 
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4 c 4 r 



3r z 



Gm 



Gm 
77 + 4 r ;>r 
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r 



(2 + 297/) h47/r 



Gm ( 
r D = r A + |5777j 2 + 2 (1 + I277) 
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8c 1 
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4 c 4 r 
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5t/ - 2?j 2 r 2 - 3r 4 



f 2 _ 



2(1 + I777) v 2 - (4 + 3877) r 2 



Gm 



2c 4 r 



■r 2 



15 (t; 2 -r 2 ) 77 + (1 + 277)^| 



(4.6a) 
(4.6b) 

(4.6c) 
(4.6d) 
(4.6e) 

(4.6f) 
(4.6g) 



The subscript 'D' denotes quantities in the De-Donder ( harmonic ) coordinates. Note 
that in all the above equations the differences between the two gauges are of the 2PN order. 
As there is no difference between the harmonic and the ADM coordinates to 1PN accuracy, 
in Eqs. ( |4.6| ), for the 2PN terms no suffix is used. The 2PN extension of the evolution of the 
orbital elements thus requires more technical care than the 1PN case due to the differences in 
the ADM and harmonic coordinates given by Eqs. ( |4.6|) . Finally using the above equations 
we have verified that the expressions given by Eqs. ( |2.8|) , relating the individual locations 
of the two bodies to the centre of mass coordinate are consistent with the corresponding 



choice in ADM coordinates, given by Eqs. (A5) - (A8) of |p6 |. 
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C. 2PN corrections to < d£ /dt > and < dj/dt > 



Starting from Eqs. fl3.5| ) and ( ^.9| ) for the far- zone fluxes in the harmonic coordinates we 
use Eqs. ( |4.6[ ) , to obtain dS / dt and djj dt in the ADM coordinates. For economy of presen- 
tation, we write the results in the following manner, (FIux)a = (Flux)o+ 'Corrections 1 , 
where (Flux) a represent the far- zone flux in the ADM coordinates. (Flux)o is a short 
hand notation for expressions on the r.h.s of Eqs. G3.5|) and ( p.9|) , where v 2 ,r,r are the 
ADM variables v A ,rA,rA respectively. For example, the Newtonian part of (dS/dt)o will 
{12v\ — llr A }. The 'Corrections' represent the differences at the 2PN order, 



be 



G 3 



15 c 5 



that arise due to the change of the coordinate system, given by Eqs. ( f4.6|) . As the two 
coordinates are different at the 2PN order, the 'Corrections' come only from the leading 
Newtonian terms in Eqs. ( |3.5|) and ( |3.9|) . 

Cm 



d£\ /dS\ G*m 3 fi 2 
dt J a \ dt J o 



15c 9 r| 

360v A - 18A0v 2 A r 2 A + 1424r A 

fdj_\ fdj\ G 3 mV(L j v) A f 
\dt) A \dt)o 5c 9 r A I 

+(363v A f A - 50^ - 363r A ) n 



(48 + 336n)v 2 A - (36 + 232n)r 2 A 
V 



Gtyi 

(4 + 68??X - (8 + 7Qn) + (2 + 82n)r 2 A 



(4.7a) 

Gm 

ta 

(4.7b) 

Note that all the variables on the r.h.s of Eqs. ( f4.7| ) are in the ADM coordinates. In the 
circular limit energy and angular momentum fluxes are again related as in Eqs. ( |3.15| ), via 
the corresponding l v 2 ' in the ADM coordinates given by 



Gm L .„ . Gm 1 2 

U - (3 - »7 - — + -(42 - 577 + 877 2 

ta I c z ta 8 



(4.8) 



From this point onwards, in this section, we work exclusively in the ADM gauge and 
hence we drop the subscript 'A' for the ease of presentation. We now have all the ingredients 
needed to calculate the 2PN corrections in < dS/dt > and < dj /dt >. We explain in detail, 
the procedure to compute < dS/dt > and only display the final expression for < dj /dt >, 
as the procedure is the same in both the cases. Starting from Eqs. ((4.7|), (|3.5|), and ( p.9| ) 
which give the far zone fluxes as functions of v 2 , r 2 , and Gm/r, we use the 2PN accurate 
orbital representation, to rewrite dS/dt as a polynomial in (1 — e r cosm) -1 . This polynomial 
is of the form 
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a N (E, h) 



(1 -e r cosuY N+1 ) 



(4.9) 



where for the convenience we have factored out du/ndt given by 
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(1 — e r cosu) 



(4.10) 



It is a straightforward algebra to show that the coefficients on(E, h) in Eq. ( |4.9| ) take the 
form 
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Gc 5 



-E) 5 fa(E,h) 



(4.11) 



where 0n(E, h) for N = 1, 2, ... 8 are given by 
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To the 1PN order Eqs. Q4.12 ) agree with Eqs. (4.15) of ||46|| . The far-zone energy flux 
(dS/dt) is a periodic function of time with period P = 2ir/n. Averaging (dS/dt), given by 
Eqs. ( |4.9| ), ( |4.11| ) and ( |4.12| ) over one time period P, we obtain 
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The integrals in Eq. ( |4. 13| ) are the Laplace second integrals for the Legendre polynomials |63 
which yield, 
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where Pn is Legendre polynomial. Using Eq. ( |1.14j ) in Eq. ( |4.13|) , we obtain < dS/dt > in 
terms of E and e r : 
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Following exactly a similar procedure, we obtain the 2PN correction to < dj jdt >. The 
final result we obtain is: 



27 



dj 

< —r- > 
dt 



4 jin (-2E)2 

5 ^(1-e 2 ) 3 
i-E) 
168 c 2 
{-E)'< 



el 



7e 



(292 + 705 677) + (197 38 + 144 34 77) e 2 + (127 + 133 77) e\ 



(1 - e 2 



- 240 - 96r/ - (30 - 12r])e 2 r - (210 - 84r/)e; 



'1 - e?l 



299 623 220 25 



1134 



252 



T] + 



351 



/1316273 815 597 292 07 



V 864 
290113 3 



336 
124403 



V 



2 \ 2 
V I e; 



96 

7187 2 > 4 
-*7 I e„ + 



7526 2869 



143 5 



6048 96 ' 48 ' ) ~ r ' V 63 224 ' ' 96 

To the 1PN order, Eqs. ([4.15 ) and (|4.16 ) agree with |]46|J48[| as required. For the special 
case of circular orbits, e r = and we observe that, < dS/dt >= uo < dj /dt > to the 2PN 
order, where u, the mean angular frequency of the relative motion, defined by u = n(l + k) 
is given by 
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It is not very difficult to trace the origin of the two types of terms in Eqs. Q4.15Q and 
( 1.16Q at the 2PN order. It is related to the fact that 'Corrections^ in Eqs. ([4.7|) , arising 
from the transformation equations connecting the harmonic and the ADM coordinates have 
a different functional form than the 2PN contributions to the corresponding far-zone fluxes 
in the harmonic coordinates. For example, in the far-zone energy flux, 1 Corrections'' contain 
a common factor (G 4 m 3 /r 5 ), unlike the 2PN contributions in harmonic coordinates which 
have only (G 3 m 2 /r 4 ) as the common factor (c.f Eqs. (|3.5| ) and ( [4.7| )). These different 
functional forms, after the averaging procedure give rise to the two different types of terms 
in Eqs. (|47T5| ) and ( ^T6|) . 

We display below < dS/dt > and < dj /dt > in terms of Gm/a r and e r , which can 
easily be obtained from Eqs. ( |4. 15| ) and ( |4.16j ), using E written in terms of Gm/a r and e r 
to the 2PN order. The required equation for E is obtained from Eqs. (|4.3|) for a r and e r by 
inverting them for E and h 2 respectively order by order. Eliminating h 2 from the expression 
for E we finally get, 
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where ( = Gm/c 2 a r . Using the above expression for E, Eq. ( |4.15|) becomes 
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while Eq. ( |4.16|) gets transformed to, 
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We observe that in the test particle limit ( 77 — > 0) and for small radial eccentricities, 
Eqs. ( PIP and QODP become 
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Such expressions for average energy and angular momentum fluxes for a test particle moving 
in a slightly eccentric orbit around a Schwarzschild black hole have been obtained by Tagoshi 
|24| , using the black hole perturbation methods: Eqs. (4.9) and (4.12) of [|4| (with q — 0). 



They are given by 
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where v and e refer to the radial velocity and the eccentricity in Schwarzschild coordinates. 
Eqs. ( f4.21|) and ( |4.22| ) are consistent, if the ADM variables a r and e r are related to the 
Schwarzschild variables v and e by 
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As stressed by Tagoshi, the fluxes reveal the more familiar coefficients in terms of a parameter 
v', related to the angular frequency in the coordinates rather than v , which is adapted to 
the radial coordinate r. For slightly eccentric orbits, v and v' are related by 
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In terms of v' the far- zone fluxes for a test particle in Schwarzschild geometry, Eqs.( |4.22| ) 
may be written as 
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In this form at the Newtonian order, one recovers the results of Peters and Mathews [f42 . 



The quantities a r and e r in the ADM coordinates are related to v' and e by the following 
relations 
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The above relations may be rewritten, in terms of the conserved energy E using 
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which are the generalizations of similar 1PN relations in ES 
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D. The evolution of the orbital elements 



In this section, we compute the 2PN corrections to the evolution of orbital elements due 
to the emission of gravitational radiation. We describe the procedure to compute the rate of 
decrease of the orbital period of the binary in some detail and display the final expressions 
for the rate of decay of other elements namely, < da r /dt > and < de r /dt >. Employing the 
heuristic argument, based on the energy and the angular momentum conservation to the 
2PN order, the rate of decrease of the orbital period, P of the two compact objects moving, 
in quasi-elliptical orbits is computed. The 2PN accurate orbital period, P = 2ir/n given in 
34]-p6[ reads as 
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Differentiating Eq. fl4.29|) with respect to t and equating dE/dt to (— < d£ j 'dt > / ' ji) and 
dh/dt to (— < dj /dt > /Gmii) we find 
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Note that, in the above equation we need < dj /dt > to the Newtonian accuracy only. 
Using in Eq. ( |4.30| ) , < dS/dt > given by Eq. ( [4.15| ) and the Newtonian part of Eq. ( [4.1 6| ) 

for < dj /dt >, we get 
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Finally inserting the expressions for e 2 and Gm/a r in terms of E and h 2 in Eq. (|4.31|) we 
obtain 
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In the expression above, P is given as a function of the masses and of the 2PN-conserved 
energy and angular momentum. This expression for P is independent of the coordinate 
system used to derive it. Since P is a measurable quantity, one would have liked to express P 
in terms of other directly observable parameters like the orbital period and some convenient 



eccentricity as in the 1PN case [|46 |. However at present, to 2PN accuracy we do not have 
any such suitable and convenient choice and therefore we leave the expression for P in terms 
of the 2PN accurate E and h 2 . 

Similarly, using the definition of a r and e r in terms of E and h 2 and following the method 
described above, we obtain after a rather long but straightforward calculation 
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To 1PN accuracy we recover the results of [13]. For the special case of circular orbits 
< da r /dt > takes the simple form 
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Eq. ( |4.35[ ) is consistent with the expression for r given in in [5E\, after taking due account 
of the coordinate transformations required to relate the ADM and the harmonic gauges for 
the circular orbits. 



V. THE 2PN CONTRIBUTION TO THE WAVEFORM 

In this section, we compute the instantaneous part of the 2PN accurate gravitational 
waveform i.e., the transverse - traceless (TT) part of the 2PN accurate far-zone field for two 
compact objects of arbitrary mass ratio, moving in a general orbit. It is given by jl6fl : 
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where R is the Cartesian observer-source distance and iV a 's are the components of N = X/i2, 
the unit normal in the direction of the vector X, pointing from the source to the observer. 
The transverse traceless projection operator projecting orthogonal to X, is given by 
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Evaluating the appropriate time derivatives of the multipole moments and performing the 
relevant contractions with N as required by Eq. ( f5.1|) , some details of which are given in 
Appendix A, we obtain 
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-(1242 -417077+ 193077 



2 Gm 



n (i v i) 



+ 



1230 - 38 1 07/ - 90t/ z 



rv 



1 



+ 60< N "» 



, Gm 



117 - 4807/ + 540t7 2 )t; 4 - (630 - 28507/ + 40507/ 2 )7; 2 f 2 



G TTt 

-(125 - 74077 + 9OO77 2 ) 7j 2 



+ (105 - 10507/ + 315077 2 )f 4 + (2715 - 85 807/ + 1260t/ 2 ) 



Gm 



-(1048 - 31207/ + 240t/ 2 ) 



G 2 m 2 



+ 



(216 - 13807/ + 4320t7 2 )t/ + (1260 - 3300t/ - 36007/ 2 )r 2 



(39 52 - 128 607/ + 36 607/ 



2 ,Gm 



r nuv 



d v j) 



12 - I8O77 + 11607/ 2 )t; 2 + (1260 - 3840t/ - 7807/ 2 )r 2 



-(664 - 23607/ + 1700t/ 5 



1 

60 



Gm 
) 

r 

2\„.4 



(66 - 15t/ - 125t/> 



Gm 



+ (90 - I8O77 - 4807/ 2 )7/r 2 - (389 + 1030t/ - 110t/ 2 ) tj~ 

r 
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+ (45 - 225?7 + 225r/ 2 )r 4 + (915 - 14407/ + 720t7 2 ) 



G m 



+ (1284 + 1090ry) 



G 2 m 2 



Gin 



Tlij 



132 + 540r/ - 580t]' 2 )v z + (300 - 1140/7 + 30077 2 )r 2 



+ (856 + 4007/ + 70077 2 



G m 



Gm 



■ r nav 



(45 - 31577 + 585t7> 4 + (354 - 21077 - 550t7 2 



Gm 



-(270 - 3077 + 27077 2 

-(638 + 140077 - 13077 2 

The "tail" contribution reads 
2G 2Gm 



Gm 



G 2 m 2 



(^fcm)tail 



C 4 R C 3 



-v„ 



ijkm 



T 



3c 

where we have used for simplicity the notation 

b^be- 11 ' 12 , b 2 = be- 97 ' m 



+ YMw,) £ ^ NbJ ^ {TR - T 



( 4 ), 



b-, = be~ 7/6 . 



(5.4e) 



(5.5) 



(5.6) 



We do not discuss the "tail" terms in this paper. Some details of these tail terms may be 



found in f|l6|,|l8[ . 



The first check on the above waveform is its circular limit, which matches with the 
wave form computed earlier in [T(J. The next check of the waveform in the general case is 
performed by computing the far-zone energy flux using 

dS c 3 R 2 



uTT h TT 

km km 



dfl(N). 



(5.7) 



dt 32%G 

The expression for dS/dt thus obtained is identical to the far- zone energy flux directly ob- 
tained from multipole moments Eq. ( |3.5| ). Of course, these checks do not uniquely fix the 
expressions in Eq. ( |5.4| ) and equivalent expressions are possible leading to the same trans- 
verse traceless parts as discussed below. 
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The above expressions for the waveform, computed using STF multipole moments differ 
from the corresponding expressions obtained by Will and Wiseman (Eqs. (6.10), (6.11) of 
||18||), using the Epstein- Wagoner multipole moments at 1.5PN and 2PN orders. Though 
the two expressions are totally different looking at these orders, even in the circular limit, 
it is possible to show that they are equivalent. The equivalence is established by showing 
that the difference between the two expressions, at 1.5PN and 2PN orders has a vanishing 
transverse-traceless, (TT) part. The easiest way of verifying this is to show that the 'plus' 
and 'cross' polarizations of the difference in the two expressions vanish at 1.5PN and 2PN 
orders f64|. In appendix B, we present the difference - at 1.5PN and 2PN orders -, between 



our waveform expression computed directly using the STF multipoles and the Will- Wiseman 
one computed using the EW multipoles and verify their equivalence. Finally we note that 



the statement in the appendix E of [pl| should more precisely read that, STF multipole 
moments presented there yield an expression for the waveform equivalent to Eqs. (6.10) and 
(6.11) of ||18|| , and not identical to it |64 |. 



VI. CONCLUSION 

In this paper using the BDI approach, we have computed the 2PN contributions to the 
mass quadrupole moment for two compact objects of arbitrary mass ratio moving in general 
orbits. Using this moment we have computed the 2PN contributions to the gravitational 
waveform and the associated energy and angular momentum fluxes. These expressions have 
already proved useful in the computation of the 2PN radiation reaction, i.e the 4.5PN terms 



in the equations of motion | 55| , using the refined balance method proposed by Iyer and Will 



56| , |57[| . Work is in progress |65| to obtain the higher order corrections to the far- zone linear 



momentum flux from the gravitational waveform presented here, extending the treatment 
of Wiseman ||66|| . It should be noted that 2PN corrections to the linear momentum flux can 
be computed only if one knows hj£ to 2.5PN accuracy. Using the 2PN accurate generalized 
quasi-Keplerian representation for elliptic orbits, we have computed here the instantaneous 
2PN contributions to < dS/dt > and < dj/dt >, the fluxes averaged over one orbital 
timescale. This is used to compute the evolution of the orbital elements, in particular P, 
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e r and d r . The method employed to compute < d£/dt > and < dj jdt > could also be 
adapted to the case of hyperbolic orbits to generalize the work of Simone, Poisson and Will 
on the head-on collision pD| . 
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APPENDIX A: STF TENSORS AND FORMULAS FOR THE WAVEFORM 

COMPUTATIONS 

We present details of the scheme, employed to compute the contributions to hjk from 
various multipole moments, as required by Eqs. ( |5.1[ ), ( |5.3| ) and ( |5.4| ). Our scheme proceeds 



in steps. In the first step, we write down schematically, the form of the desired time derivative 
of the STF multipole moment, using the compact notation {}, introduced by Blanchet and 
Damour |TJ]. Here {} denotes unnormalised minimum number of terms, required to make 
the expression symmetric in all the indicated indices. The second step involves peeling, where 
by observation and counting, we rewrite the expression obtained in the step 1, as STF on 
the free indices - i and j in our case -. In Step 3, we contract, the final expression of step 



2 with appropriate number of N's as required by Eq. Q5.1|) . The actual evaluation of the 



result of step 3 is performed on Maple [|53|| . In all the formulae, Sl, denotes the symmetric 

r^j if the object is and Sl = 



version of the object under consideration; e.g. Sl = if the object is and Sl = j\ m: 



if the object is J L ; - the object in the formula is obvious from the context. 
The unnormalized symmetric blocks . 

${ijSa} = $ij S a + 5i a Sj + 5j a Si , (Ala) 
U{ijS a b} — fiijSab + ^iaSjb + 5ibS a j + 

8jaSib + SjbS a i + 5 a bSij , (Alb) 
5{ij5 ab } = Sijdab + 5 ia 5 jb + 5 ib 5 aj , (Ale) 
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STFjj a (ijj ) = Sij a — —S^ijS a -ju , 

STFij ab (Iij ab ) = Sij a b — -S{ijS a b}tt + 7^:fi{ij8ab}S S stt , 

STFjj a fr c (/jj a 5 c ) Sijabc — 8{ i j S a bc}tt ~\~ ~^7^{ij$abSc}sstt j 

1 1 
^^-F ij abcdi^I ij abed) Sij a b c d {ij ^ abcd}pp ~^^${ij$abS c d}ttqq 

~~^^fi{ijfiab3cd}Sppqqtt ■ 



(A2a) 
(A2b) 
(A2c) 

(A2d) 



The 'Peeling'. 

STFjj a (/jj a ) = STFjj l^jo — g^m'Sjtt| ) 



STFijabil'ijab) — STFjj-j Sij a b — — 



^iaSjbtt + ^ibSjatt + ^baSjitt 



+ 



2 

35 



$ia$jbSttss 



(A3a) 



(A3b) 



STFjj a 5 c (/jj a ^ c ) STFjj ^SV/afic g 2 ^SiaSjbcpp $ibSjacpp ~\~ ^icSabjpp^j 
~t~ i^abSijcpp ~\~ $ acS ijbpp ~\~ $bcSijapp^J 



"I - ( &ib^jc$appqq ^ia^jc^bppqq ^ia^jb^eppqq 

l 

IT 



(A3c) 



1 

+ 99 



^ijabed YY ^ ( ^i a ^jbcdpp $ibSjacdpp $icSabjdpp &id$abcj pp 
~t~ {^adSbcijpp $abSdcijpp 

~^O ac Sbdijpp $bcS adijpp ^bd^acijpp ^cd^baijpp 

2 ia$ jbS cdppqq $ia& jc^bdppqq 
$ic$ jbS adppqq ^ia^jdScbppqq 
$ id$ jbS acppqq ~\~ $ic$jdSbappqq 

+2^S ic 5 a b + 5 ia 5 c b + 5ib5 ca jSjdppqq + 2[5id5 a b + 5 ia 5db + 5ib5dajSj cppqq 
+2 (dialed + 8ic8da + Sid^ca^ Sjbppqq + 2 (bibbed + ^ic^db + S id 5 b ^j S. 
+ I ScdSab + 5 ca 5 d b + ^da^cft ) Sjippqq 



J jappqq 



693 



ftiabjbbcd + S ia Sj c S b d 
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ac ) &jd 

The contractions with N 

(3) f (3) (3) 1 

STF ija (I ija ) N a = STFijtS ija N a - - NiS jtt > , 



S 



ppqqtt 



STF ijab{l\fab) N ab — STF^ ^S\f ab N ab - - 
r(5) 



r 6 l 

STF ija6c (/^ 6c )A abc = STFj,|s^ 6c A abc - -A^S 1 ]^ - - 



+ 7^NijSt t g S 



ijtt 

- c ? ( ' 5 " ) A/" 4- 
2 '-'ijcpp 1 v c ^ 



ic( 5 ) TV- 4- —TV - 9^ 1 

63 63 ya a-PPIQ J ' 

'ijabcd(Iijibcd)Nabcd = STF ij^S^ltejNabcd — JjA, 



?(6) 



j ^ J y «bcd ^jbcdpp ^ ^ ^ijcdpp Iy cd 



693 ^ppqqttj ■ 
The current multipole moments. 

£pq(i'Jj)pL STFjj | e p gjJ, pL | 

e pq{iJj)p Nq = $TFij^€pqiSj p Agj- , 



+ 1^AT ^ {6) 4-— /V <?( 6 ) 4-A^ 6 ) 

99 ^ c< ^ cdppqq 99 i d ppqq 99 vppqq 
_il/V-.9 {6) 1 



e pq(ijj)l a Nqa = STF je pg i 
e pq(iJj)labNqab = STFjj|e pg j 
e pq(iJ j)pabc N qabc = STFjj<e p9 j 



(3) 



jpa - 1 * ?a 



ptt Jv gj 



1 



J jpab ly qab ~Z I - 'A - ' . 



(5) 



Sjpabc Nqabc g ( 'S'pfcctt ^qjbc + SjpcH - N </- 



<7 (4) AA 

(5) N 



1 
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_i 1_ I c( 5 ) /\r . 

l ^pttvv 1 ' qj 



(A3d) 



(A4a) 
(A4b) 

(A4c) 



(A4d) 



(A5a) 
(A5b) 
(A5c) 

(A5d) 



The explicit computations of the above equations require the following identities, which are 
easily derived, using the rules governing the product of e's. The identities are 



STF,,- 


£pqi NqVjLp^ 


= STFj, 


STF;, • 


£pqi NgVjLp^ 


= STF,,- 


STFj, ^6pgj AAqjZ/pj 


= STFj, 


STF,,- 


^pqi NqUpLj J 


= STF;.,- 


STFy 


£pqi NgVpLj J 


= STFy 



-(N.v) !/lj + (N.n)rp i j, (A6a) 

-(N.v)^ + (N.n)r % J, (A6b) 

(N.n) r NiVj - (N.v) y^j , (A6c) 

-(N.v) yij + (N.n) rytVj + (rr)N iyj - r 2 A^j , (A6d) 

(N.n) r Vij - {N.v)yiVj - rr N^j + v 2 A^y, j , (A6e) 
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STF J e w iV> p (L.N) 



STFij \e pqi N qj y p (L.N) 



STFiA e pqi N q y jp (L.N) 



STFjJ e pqi N q yjV p (L.N) 



STFjJ e p<?i N q y p Vj(L.N] 



STFiA e pqi N q v jp (L.N] 



where L p = e pk iy k vi. 



STF^j^ 2 - (N.v) 2 ]A i% - + [(N.n)(N.v) - r]riV^ 
+ [r(N.v) -i; 2 (N.n)]r^ 
STFijjfr - (N.n)(N.v)]rJV i% + [(N.n) 2 - l]r 2 N iVj + 
[(N.v) - r(N.n)]r 2 iV, 



STF fi |[r - (N.n)(N.v)]r^ + [(N.n) 2 - \}r 2 y iVj + 
[(N.v) - r(N.n)}r 2 N iyj 



STFiA [v 2 - (N.v) 2 ]^- + [(N.n)(N.v) - r\r Vl v d 



+ [r(N.v) - v 2 (N.n)]riVj^ 
STF -|[r - (N.n) (N.v)] r + [(N.n) 2 - l]r 2 w 
+ [(N.v) -r(N.n)]r 2 iV> 



j r ' 



STF y |[v 2 - (N.v) 2 ]?/^ + [(N.n) (N.v) - r}rv 
+ [r(N.v) - v 2 (N.n)]rN i v j \ , 



(A6f) 



(A6g) 



(A6h) 



(A6i) 



(A6j) 



(A6k) 



APPENDIX B: THE EQUIVALENCE TO WILL-WISEMAN WAVEFORM 



The expression for the gravitational waveform, obtained by Will and Wiseman |1J| differs 
from our waveform expression at the 1.5PN and the 2PN orders. We give below the difference 
in the waveform expressions at these orders and show that the two polarization states, h + 
and h x of the difference are zero at 1.5PN and 2PN orders. 

-(N.v) (N.n) 2 Vij+Qrn^Vj) 
+ (N.n)(N.v) 2 2 n {i v i} + 3 f n„ 



{(H 



TT\ (1-5) (h TT\(l-5) 
krnJBDI {"'km) WW J 



-(N.v) 3 ^- + 3 (N.n)r 



v 2 nij + Vij - 2rn {i Vj) 



(N.v) 



v 2 nij + - 2rn {i Vj) 



(Bla) 
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Uh TT ) {2) 



(h 



TT\(2) ! 
km) WW J 



1 _ Gm I 

Vijkm — <, (1 - 5?7 + 5?7 



15 c 



r 

i I / 9„.2 



12(N.v) 4 n 



-3 (N.n) 4 [3v 2 - 15r 2 + 



+6(N.n) 3 (N.v) 



3w 2 - 15r 2 



-6(N.n)(N.v) 3 (9rn ii + 4n {i v 



-3(N.n) 2 (N.v) 2 
-(N.v) 2 



3v 2 - 15r 2 



Gm 



ru 



36fri(iVj) 



-(39 - 125?7 - 5?7 



(51 - 185r? + 55r/ 2 )t> 2 - (117 - 375r? - 15r/ 2 )f 2 
Gm s 



-24 I 1 - 5r] + hrj 2 )rn {i v j) + 12 ( 1 - 5r? + hrf )v i:j 



+2 (N.v) (N.n) 



27 ( 1 — 5?7 + hrf ) f t> 2 n i3 



+ (39 - 125r/-5r/ 2 )[w 2 - 



Gm, 



-(171 - 645r7 + 255ri 2 )f 2 jn^Vj) + 27 ( 1 - 5r? + hi] 1 ) f u, 

-(N.n) 2 

2 Gm 
— 3t> 



1 - 5?7 + 5r? 2 ) ^-9 + 45rV 
(riij - 2rn{iV j )) 



+ ^(30 - 80?7 - 50 r] 2 )v 2 - (72 - 150?? - 240 rf)r 2 



-(42 - 140r7 + 10?7 2 



Gm 



(39 - 125r/ — 5r/ 2 ) ( 
(v 2 riij - 2 r n^Vj) + 



Gm 



v 2 ) + (117 - 37577- 15 r/ 2 )r 2 



(Bib) 



The two independent polarization states of the gravitational wave h + and h x are given 
by h + = \{piPj — <li<lj \ hJT and h x = \{pi<lj + Pj q^j hj? , where p and q are the two 
polarization vectors, forming along with N an orthogonal triad [ IB jj7 , n| . Note that there 
is no need to apply the TT projection before contracting on p and q. Consequently, we 
write the difference in the waveform at the 1.5PN and the 2PN orders as 



{{hj^ww ~ (/'// )/;/)/! = CiVij + (211^ + C^niiVj) 



(B2) 



The polarization states h + and h x , for Eqs. 



are given by, 
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K = - hi Pj - ft qj ) ( Ci Vij + C2 riij + C3 n^Vj) 



|((p.v) 2 - (q.v) 2 ) + § ((P.n) 2 - (q.n) 2 ) + | ((p.n) (p.v) - (q.n) (q.v)) , (B3a) 
2 (ft 1j + Pj ft) (Ci + C2 riij + C3 

d (p.v) (q.v) + C 2 (p.n) (q.n) + | ((p.n) (q.v) + (p.v) (q.n)) (B3b) 



For the explicit computation of Eqs. (p3|), we use the standard convention adopted in 

1161 , 1671 , 1181 , which gives, p = (0, 1, 0), q = (— cosi, 0, sini), N = (sinz, 0, cosi), 

n = (cos (j), sin 0,0), and v = (r cos <p — ru sin 0, r sin <p + r u cos 0,0), 

where n and v are the unit separation vector, and the velocity vector respectively, is the 

orbital phase angle, such that the orbital angular velocity u> = d(f>/dt and £ i' is the inclination 

angle of the source. 

A straightforward but lengthy computation shows that h + and h x , given by Eqs. ( |B3| ) 
vanish, both at the 1.5PN and the 2PN orders. This establishes the equivalence of our 
waveform expression, Eqs. (5.3) and (|5.4Q with the WW one given by Eqs. (6.10) and (6.11) 
of PH. 
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